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ABSTRACT 

In  this  paper  we  develop  a  comprehensive  framework  for  the  study 
of  decentralized  estimation  problems.  This  approach  imbeds  a  decen¬ 
tralized  estimation  problem  into  an  equivalent  scattering  problem, 
and  makes  use  of  the  superposition  principle  to  relate  local  and  central¬ 
ized  estimates.  Some  decentralized  filtering  and  smoothing  algorithms 
are  obtained  for  a  simple  estimation  structure  consisting  of  a  central 
processor  and  of  two  local  processors.  The  case  when  the  local 
processors  exchange  some  information  is  considered,  as  well  as  the  case 
when  the  local  state-space  models  differ  from  the  central  one. 
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I.  Introduction 

The  study  of  coop lex,  large-scale  systems  involves  frequently  the 
solution  of  decentralized  estimation  problems .  This  is  the  case  for 
example  when  one  considers  some  systems  which  are  inherently  decentral¬ 
ized,  such  as  power  systems  [1] ,  or  when  communications  between  sub¬ 
systems  are  limited  or  unreliable,  as  in  military  command  and  control 
systems  [2].  Also,  in  seme  applications,  the  volume  of  information  to 
be  processed  is  so  large  that  a  decentralized  or  parallel  processing 
structure  is  required.  Such  structures  have  been  proposed  for  example  in 
image  processing  (if  the  strip  Kalman  filter  discussed  by  Hoods  and 
Radewan  [3] )  or  in  order  to  estimate  the  location  of  moving  sources  in  a 
distributed  sensor  network  [4] . 

An  additional  incentive  to  consider  decentralized  estimation  struc¬ 
tures  is  the  recent  progress  in  VLSI  (very  large  scale  integration) 
technology  which  makes  it  now  feasible  to  implement  estimation  algorithms 
in  a  parallel  fashion  IS].  However,  in  spite  of  these  potential  applica¬ 
tions,  it  is  only  recently  that  decentralized  estimation  problems  have 
been  the  object  of  seme  attention.  This  is  due  in  part  to  the  complexity 
of  these  problems:  to  describe  entirely  a  decentralized  estimation  struc¬ 
ture  consisting  of  interconnected  local  processors  one  needs  to  specify 
for  each  processor  a  local  model,  a  set  of  measurements  and  a  local 
objective  function  (see  Barta  {6]  for  a  general  formulation  of  decentral¬ 
ized  estimation  problems ) .  In  addition,  one  needs  to  describe  the 
topology  and  the  capacity  of  the  cossnunications  network  connecting  all 
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processors.  Thass  various  elements  aim  at  capturing  the  fact  that  in 
a  decentralized  estimation  structure,  each  processor  by  itself  has  only 
a  limited  knowledge  of  the  physical  process  under  study,  so  that  to 
overcome  this  limitation,  the  local  processors  must  communicate.  Also, 
even  if  all  processors  cooperate,  their  objectives  may  not  be  the  same. 

Another  important  feature  of  decentralized  estimation  structures  is 
that  they  provide  usually  more  than  one  way  to  solve  a  given  problem. 

This  explains  why  several  approaches  have  been  proposed  to  study  decen¬ 
tralized  estimation  problems.  Barta  [6]  has  used  a  game- theoretic  ap¬ 
proach  that  allows  for  a  large  amount  of  flexibility  in  the  choice  of 
information  flows  between  local  processors.  By  contrast,  Sanders  et.al. 

[7]  and  Tacker  and  Sanders  [8]  have  considered  decentralized  estimation 
filters  with  a  fixed  structure.  Another  approach,  which  was  followed  by 
Speyer  19],  Chong  [10]  and  Willsky  et.al.  [11]  is  based  on  decomposing  a 
central  estimation  problem  into  smaller,  local  ones.  This  approach  will 
be  at  the  center  of  our  discussion  here. 

The  objective  of  this  paper  is  to  propose  a  more  systematic  framework 
for  the  study  of  decentralized  estimation  problems.  To  do  so,  we  will  use 
the  scattering  framework  for  linear  estimation  introduced  by  Kailath  and 
his  coworkers  [12]- [15].  In  this  context,  to  every  local  or  centralized 
linear  estimation  problem,  we  can  associate  an  equivalent  scattering  problem 
which  is  specified  by  the  knowledge  of  soste  scattering  parameters  and  of  a 
set  of  internal  sources  for  the  scattering  medium.  Then,  if  the  scattering 
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parameters  describing  the  various  local  or  central  estimation  problems 
can  be  transformed  so  that  they  give  rise  to  the  same  media,  we  can  use 
the  superposition  principle  for  linear  scattering  media  to  relate  the 
local  and  central  estimates.  This  technique  will  be  used  here  to  obtain 
several  decentralized  estimation  algorithms.  To  illustrate  these  results 
we  will  also  consider  several  random  fields  problems  of  the  type 
discussed  in  [11] . 

This  paper  is  organized  as  follows.  In  Section  II  the  scattering 
framework  for  linear  estimation  problems  is  introduced,  and  in  Section  III 
a  brief  review  of  the  properties  of  scattering  media  is  given.  Then,  in 
Section  IV  we  obtain  some  decentralized  smoothing  and  filtering  algorithms 
for  a  simple  decentralized  structure  consisting  of  a  central  processor  and 
of  two  local  processors  which  do  not  communicate  between  each  other.  When 
these  local  processors  are  allowed  to  communicate,  the  results  of  Section  IV 
need  to  be  modified  as  shown  in  Section  V  where  we  consider  several  problems, 
the  smoothing  update  and  the  real-time  smoothing  problem,  which  appear  in 
this  context.  The  case  when  the  state- space  models  available  to  the  local 
processors  differ  from  the  global  model  is  also  considered  in  Section  VI. 
Finally,  in  Section  VII,  we  discuss  some  possible  extensions  of  this  work 
to  networks  where  delays  or  noise  can  appear  in  the  communication 
channels  between  processors. 

These  constitute  only  a  few  of  the  issues  appearing  in  the  context 
of  decentralized  estimation.  In  a  companion  paper,  the  scattering 
framework  will  also  be  used  to  obtain  some  decentralized  algorithms  for 
weakly  coupled  and  singular  perturbed  systems. 
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II.  The  scattering  framework  for  filtering  and  smoothing  problems 

The  scattering  framework  for  the  study  of  linear  estimation  problems 
has  been  discussed  in  great  detail  in  112] -[15].  Me  will  therefore  give 
here  only  a  brief  presentation  of  the  main  aspects  of  this  theory.  We 
consider  the  state-space  model 

x(t)  -  Ax (t)  +Bu(t),  0<t<T  (2.1) 

with  observations 

y(t)  -  Cx (t)  +  v (t)  (2.2) 

where  u(t)  and  v(t)  are  same  white  noise  processes  uncorrelated  with 
the  initial  conditions  x(0),  i.e. 


and 


(2.3) 


(2.4) 


with  R>R'>0  (the  estimation  problem  is  assumed  to  be  nonsingular). 
The  a  priori  information  on  the  initial  conditions  is  given  by 

B[X(0)1  -  x(0),  E[x(0)x'(0)]  ■  P(0)  (2.5) 

m  A 

where  x(0)  ■  x(0)  -  x(0).  in  the  following,  it  will  be  convenient  to 
view  this  information  as  an  additional  observation  on  x(0)  given  by 
m  ■  x(0)  ■  x(0)  +  w 


(2.6) 
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where  E[w]*0  and  E(ww' ]»p(0) .  Then,  if  we  introduce  the  Hilbert 
spce 

X  -  H(x(t) ,  0<t<T),  XQ  -  H(x(0) ) 

U  -  H(u(t) ,  0<t<T) ,  V  -  H(v(t) ,  0<t<T) 

X  -  H(y(t),  0<t<T) 


the  space  y  of  observations  is  included  in  the  space  G  ■  x0@U©V  =  X@V 
which  specifies  the  total  information  contained  in  the  system  (2.1)-(2.2), 
i.e. 


X  C  G, 


(2.7) 


but  in  general  this  inclusion  is  strict,  so  that  from  the  observations 
X  one  cannot  recover  cong>letely  x(t)  and  v(t)  for  0<t<T. 

However,  as  was  observed  by  Weinert  and  Desai  [16] ,  the  information 
lost  by  X  can  be  entirely  recovered  by  constructing  the  dual  process 
2  ( • )  defined  by 


X(t)  -  -A'X(t)  +  C'R-1v(t),  X (T)=0 
*(t)  *  B'X(t)  +  u(t) 


(2.8) 


with  the  observations 

n  «  P (0)X(0)-w  .  (2.9) 

In  this  case,  if 


Z  •  H (n; z (t) ,  0<t<T) 

denotes  the  Hilbert  space  spanned  by  the  dual  process,  one  can  show 


the  following  result 
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Lemma:  Z  is  the  orthogonal  complement  of  Y  in  G,  i.e. 

(i)  Z  J-Y  and  (ii)  Z©Y  «  G 

The  proof  of  (i)  is  obtained  by  direct  verification  (see  also  [16] ) . 
To  prove  (ii)  one  needs  only  to  note  that  the  system 


with  observations 


(2.10) 


(2.11) 


is  invertible.  The  inverse  system  is  given  by 


(2.12) 


(2.13) 


with  the  boundary  conditions 

X(T)»0  x(0)  +P(0)X(0)  -  a+n  . 


(2.14) 


The  solution  of  the  two-point  boundary  value  problem  (2. 12) -(2. 14) 
enables  us  to  reconstruct  entirely  G  ■  x0®  u  ©  v  from  Y  and  Z 


so  that  (ii)  is  satisfied. 
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This  result  can  be  used  to  obtain  the  scattering  framework  for  the 
smoothing  and  filtering  problems.  By  conditioning  equations  (2.12)  and 


(2.14)  with  respect  to  Y,  and  by  using  the  fact  that  Y±Z,  one  has 


with  A(t|t)=0  and 

x(C)|t)  +  P(0)A(o|t)  *  m  -  x (0)  .  (2.16) 

Here  x(t|T)  and  A (t|T)  denote  the  linear  least-squares  smoothed 
estimates  of  x(t)  and  A(t)  given  m  and  y(s)  ,  0<s<T.  Then,  by 
discretizing  (2.15),  one  gets 


and  one  can  view  (2.17)  as  obtained  from  a  scattering  medium  with 
infinitesimal  scattering  matrix  S(t+A,t),  and  interval  sources  q(t)A 


(see  Figure  2.1). 
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By  composing  together  these  infinitesimal  layers,  we  obtain  an 
aggregate  scattering  medium  which  is  described  in  Figure  2.2  (cf.  Section 
III  for  a  description  of  the  rules  of  composition  of  scattering  layers) . 
The  incoming  waves  for  this  medium  are  given  by  the  boundary  conditions 
(2.16),  and  the  solution  of  the  scattering  problem  requires  the  com¬ 
putation  of  the  internal  variables  (x(t|T),  X(t|T)).  This  can  be  done  in 
a  variety  of  ways  which  are  discussed  in  detail  in  [14],  [15].  Among 
existing  methods,  a  method  which  is  particularly  simple  is  the  Mayne- 
Fraser  two-filter  formula  [17],  [18] 

P^CtJxft  T)  =  Pp1(t)xp(t)  +  p”1  <t)acB (t)  (2.18) 

where  P  1(t)  -  P_1(t)  +  P_1(t)  and  where  (x  ,  P  )  and  (x_,  P  )  denote 
S  F  B  F  F  B  B 

respectively  the  forward  and  backward  linear  least-squares  filtering 
estimates  and  error  covariances  of  x(t).  They  satisfy  the  equations 


x  (t)  =  (A-P  C ' R~*C) X  ( t )  +  P  C'R~  y(t) 

r  F  F  i 

xp(0)  *  x(0) 

and 

-Xgft)  *  (-A-PBC'R"1C)xB(t)  +  PBC'R-1y(t) 

p"1(T)x_(T)  *  0 
B  B 

with 

•  -1 

Pp  -  APp  +  PpA'  +  BB'  -  PfC'R  CPp  ,  Pp(0)  m  P(0) 

-PB  =  -1PB  -  PBV  *  BB’  -  PsCR-1CBb  .  p;1™  -  0  . 


(2.19) 


(2.20) 


(2.21) 


x> 


A  detailed  discussion  of  this  formula  can  be  found  in  Wall,  Willsky 
and  Sandell  [19] . 

This  shows  how  to  compute  the  smoothed  estimate  x (t J T) .  To  compute 
the  filtered  estimate  x(tjt)  -  xp(t),  one  needs  only  to  use  (2.19).  By 
discretizing  (2.19)  and  by  dropping  subscripts,  one  gets 

x(t+A,t+A)  -  (I+(A-PC,R“1C)A)i(tjt)  +  PC*R~^y (t)A  (2.22) 

which  can  be  viewed  as  obtained  from  a  transmission  medium  with 

infinitesimal  transmission  matrix  T(t+A,t)»I  +  (A-PC'R  ^C)A  and 

with  internal  sources  PC'R  (t) A  (see  Figure  2.3).  These  infinitesimal. 

A 

layers  generate  an  aggregate  transmission  medium  with  initial  x(0),  and 
the  solution  of  the  transmission  problem  requires  the  computation  of 
the  internal  variables  x(t|t). 

Instead  of  solving  this  problem,  it  is  sometimes  more  convenient  to 
consider  the  information  form  of  the  Kalman  filter,  i.e. 

d(t)  -  (-A-P-1BB,)d(t)  +  C'R~^y (t)  (2.23) 

A  -1  ^  -i 

where  d(t)  »  P  (t)x(t  t)  and  d(0)  *  P  (0)x(0).  The  transmission  problem 
associated  to  this  form  has  the  advantage  that  its  infinitesimal  sources 
are  C'R  ^y (t) A,  so  that  they  depend  only  on  the  original  model  instead  of 
introducing  P(t)  as  well.  This  property  will  be  exploited  in  Section  IV 
to  obtain  some  decentralized  filtering  algorithms. 
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Figure  2.3.  The  infinitesimal  transmission  layers  associated 
to  the  filtering  problem. 
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111.  Some  facts  from  scattering  theory 

The  algebra  of  scattering  matrices  was  first  introduced  by 
Redheffer  [20] ,  [21] ,  and  a  discussion  of  the  main  aspects  of  this 
theory  can  be  found  in  [12] -[15].  Only  a  few  basic  facts  from  this 
theory  will  be  required  here. 

A  scattering  medium  is  specified  by  a  scattering  matrix 


c :) 


(3.1) 


and  by  same  internal  sources 


(3.2) 


(see  Figure  3.1).  The  parameters  a  and  b  (which  can  be  operators,  or 
matrices)  are  called  the  transmission  coefficients  of  the  medium.  They 
describe  the  portion  of  a  wave  1. (resp.  r_)  travelling  from  left  to 
right  (resp.  from  right  to  left)  which  is  transmitted  by  the  medium. 
Similarly,  the  parameters  c  and  b  are  called  the  reflection  coefficients 
and  describe  the  parts  of  the  waves  1_  and  r_  which  are  reflected.  Then, 
the  outgoing  waves  are  given  in  function  of  the  incoming  waves  by 


^  ^  ^  ^ 


(3.3) 


It  is  assumed  here  that  the  scattering  medium  is  linear. 
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to  the  filtering  problem. 
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To  compose  two  scattering  layers ,  one  needs  to  introduce  a  special 
product  different  froa  the  ordinary  matrix  product  which  is  called  the 
star  product  (see  [20} , [21] ) .  For  the  scattering  layers  described  in 
Figure  3.2,  one  gets 


(Ad-bC)*^ 
c+dC(I-bC)_1a 


B+AbCl-Cb)~1D 

d{I-Cb)-1D 


) 


(3.4) 


where  the  transmission  and  reflection  coefficients  of  S  can  be  obtain 
by  usin?  standard  flow  graph  rules  [22].  To  compose  internal  sources, 
one  needs  also  to  use  a  special  sum  called  the  assembly  sum.  For  the 
internal  sources 


of  Figure  3.2,  it  is  denoted  as  o  where 


A(I-bC)-1(q++bk_)' 
d(l-Cb)"1()crCq+) 


(3.5) 


(3.6) 


transmitted  wave 


Composition  of  two  scattering  layers. 


-17- 


Esample :  Zf  we  consider  the  scattering  medium  generated  by  the  smoothing 
problem  of  Section  II,  i.e. 


sx  -  S(t,0) 


(*<t,0) 
w(t, 0) 


-P(t,0) 

¥(t,0) 


(3.7) 


Qx  -  q(t,0)  - 


(3.8) 


the  previous  composition  rules  can  be  used  to  obtain  a  set  of  dif¬ 
ferential  equations  for  the  entries  of  S(t,0)  and  q(t,0).  Thus,  if 


I-t-AA  -BB' A 

S  -  S(t+A.t)  -  ' 

\c*B  CA  I+A’A 


and 


(3.9) 


q2  -  q(t)A  . 

\-C’R  y(t)A; 


(3.10) 


one  has 


S2  *  S1  "  S^+A.O) 


$+(A-KC)$A 
C'lf  1C0A 


P+ ( AP+PA ' +BB ' -PC ' R" XCP ) A 


mtA-KC)  'A 


U  0(A)  (3.11) 


where  K{t)  -  P(t,0)C'R  1  and  where  the  time  arguments  have  been  omitted. 


This  shows  that 
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♦(t,0)  -  (A-KC)<* 

P(t,0)  •  AP+PA,^«»,-PC,*"ICP 

»(t,0)  -  TC'r"1^ 
f (t,0)  -  f (A-KC) ' 


,  P(0,0)-0 

,  11(0, 0)-0 
,  f(0,0)-X 


(3.12) 


so  chat  P(t,0)  is  tha  solution  of  cha  Riccati  aquation  aasociatad  to 
tha  fil taring  problsa  with  zaro  initial  conditions,  # (t,0)  -  f' (t,0) 
is  tha  transition  matrix  of  tha  corrasponding  filtar  and 


w  . 

h(t.O)  -  C  (a,0)C*lf\*(a,0)ds 

JQ 


(3.13) 


is  tha  aasociatad  obsarwahility  matrix.  By  composing  aourcas,  ona 

gata  also 


Ql  °  Q2  -  q(t+0,0) 

.  I q+(t,0)  +  (  (A-KC)q+(t,0)-HCy(t)  )A 
yq_(t,0)-fC,R'1(y(t)-Cq+(t,0))A 

so  that 

q+(t,0)  •  (A-KC)q+  ♦  Ky  ,  q+(0,0)-0 

.  (3.15) 

q_(t,0)  •  -#*C'R  (y-cq^)  ,  q_(0,0)»0. 

This  shows  that 

q+(t,0)  •  xQ(t|t) 


) 


(3.14) 


(3.16) 
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where  x^(tjt)  denotes  the  Kalman  filtering  estimate  obtained  for  zero 
initial  conditions  (x(0)-0  and  P(0)-0)  and 

t 

q_(t,0)  -  -  f  4)' (3,0)C,R'*1V0(s)ds  (3.17) 

0 

where  vQ(t)  *  y(t)  -  CxQ(t|t)  is  the  innovations  process  for  zero 
initial  conditions. 


One  drawback  of  the  star  product  for  composing  scattering  layers  is 
that  this  product  is  nonlinear.  This  leads  for  example  to  some  nonlinear 
differential  equations  such  as  (3.12)  for  the  entries  of  s(t,0).  To 
linearize  these  equations,  Redheffer  noted  in  [20]  that  one  needs  only 
to  transform  the  scattering  medium  of  Figure  3.1  into  a  transmitting 
one.  This  can  be  done  by  expressing  the  waves  (r+,r_)  at  the  right  of 
the  scattering  layer  in  function  of  the  waves  (1_,1+).  By  using  flow- 
graph  inversion  rules  as  shown  in  Figure  3.3,  one  obtains 


(3.18) 


where  T  is  the  transmission  matrix  of  the  medium  and  where  R  describes 
the  transmission  sources.  Then,  the  rules  of  composition  of  trans¬ 
mission  layers  are  the  same  as  those  of  ordinary  matrix  multiplication 


so  that  if  a  medium  (T,R)  is  obtained  by  composing  two  layers 
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3.3 


Transformation  of  a  scattering  medium  into 
a  transmitting  one. 
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(Ti ,  R^ ) ,  i*l,2  one  has 

T  -  ,  R-R2  +  TJR1  . 


(3.19) 


We  note  however  that  since  scattering  problems  involve  the  specification 
of  incoming  waves  at  the  boundaries  of  the  medium,  if  we  transform  a 
scattering  problem  into  a  transmission  one,  one  gets  a  two-point  boundary 
value  problem.  Such  problems  have  been  studied  by  Reid  [23],  and  more 
recently  by  Sidhu  and  Desai  [24]  in  the  context  of  linear  estimation. 

In  the  following,  in  order  to  solve  scattering  problems,  we  will 
make  repeated  use  of  the  superposition  principle.  This  principle  applies 
only  to  linear  scattering  media.  Let  S  be  the  scattering  matrix  of  a 
linear  medium,  and  let  (r^,  l^J  be  the  outgoing  waves  of  the  medium 
associated  to  the  incoming  waves  (lSrS  and  internal  sources 


for  i*l,2.  Then,  the  superposition  principle  states  tha'.  the  outgoing 

12  2 

waves  associated  to  the  incoming  waves  (l_,r_)  =  (1+1  ,r'+r_) 
and  the  sources  Q*Q1+Q2  are  given  by 


(r+*l+) 


(r;  +  r+, 


*  £ 


(3.20) 
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IV.  The  decentralized  smoothing  and  filtering  problem 

A  simple  decentralized  estimation  problem  can  be  described  as 
follow:  let  i"l, 2  be  two  stations  which  observe  independently 
a  process 


*(t)  -  Ax(t)  +  Bult) , 


0<t«T  . 


(4.1) 


The  measureawnts  obtained  by  S .  i*l,  2  jure  denoted  as 


y.  (t)  -  C.x(t)  +  v. (t) 
i  x  x 


(4.2) 


where  and  jure  some  independent  white  Gaussian  noises  such  that 


E 


vx(t> 


L\ 


(v^(s)v^(s) ) 


5 (t-s) 


(4.3) 


Then,  if  we  assume  that  and  S2  are  linked  to  a  central  processor  P, 
and  if  we  want  to  compute  the  smoothed  or  filtered  estimate  of  x(t) 
given  ( - )  and  y2 ( - ) ,  two  strategies  are  possible: 

(i)  S,  and  send  their  measurements  to  the  central  processor 
which  computes  x,  the  estimate  of  x  given  y^(-)  and  y^t*). 

(ii)  i»l,2  processes  its  own  information  locally  and  computes 

a(1)  *(i) 

x  (the  estxmate  of  x  given  y^-)).  Then  x  i»l,2, 
is  sent  to  the  central  processor  which  in  turn  seeks  to  com- 

A (J)  A (2)  A 

bine  x  and  x  to  obtain  the  centralized  estimate  x. 

The  first  of  these  strategies  is  centralized  and  any  of  the  existing 
smoothing  and  filtering  algorithms  can  be  used  to  compute  x.  The  second 
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*(t)  -  Ax(t)  +  Bult) , 


0<t«T  . 


(4.1) 
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i  x  x 


(4.2) 
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E 


vx(t> 


L\ 


(v^(s)v^(s) ) 


5 (t-s) 


(4.3) 
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A.  Decentralized  smoothing 

In  the  following,  we  will  write  as  Y^  =  H(m,-y^(s),  0<s<T)  the 
Hilbert  space  spanned  by  the  observations  available  to  and  we  will 
denote  by 

x(t|T)  *  E lx  t)|Yi,Y2]  (4.4) 

and 

x (l)  < 1 1 T)  -E[x(t)|Y]  i-1,2  (4.5) 

the  centralized  and  local  smoothed  estimates  of  x(t). 

Then,  by  using  the  scattering  framework  of  Section  II,  the  estimates 

*  i  a (i)  > 

x(t|T)  and  x  (t]T)  can  be  obtain  by  composing  some  global  and  local 

scattering  layers  such  as  those  depicted  in  Figure  4.1a  and  4.1b.  One 

important  feature  of  these  elementary  scattering  layers  is  that  they 

correspond  to  nearly  identical  media.  The  only  difference  is  that  for 

the  global  medium,  the  reflection  matrix  associated  to  the  left  incoming 
/>  ,  -1-1 

wave  x(t|T)  is  (C'R  C  +C'R  C  )A  and  that  for  the  local  medium,  the 

X  x  X  »  *  « 

reflection  matrix  for  x ^ (tjT)  is  C^R^C^A. 

This  means  that  we  cannot  apply  directly  the  superposition 
principle  to  these  scattering  layers.  However,  a  very  simple  transfor¬ 
mation  can  be  used  to  make  the  local  medium  identical  to  the  global  one*. 

~  (1)  i 

one  needs  only  to  replace  the  local  reflection  matrix  for  x  (t|T)  by 
(C^Rj^C^+C^2 ^2 ^  ^  and  comP®nsat®  for  the  additional  reflected  wave 
C^R./CjX^  ( 1 1 T) A  by  a  fictitious  source  equal  to  -CjRj^CjX^  (t | T> A 


as  shown  in  Figure  4.1c. 
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(C1R“  C  x  (t|T)+C!R“  C  X 


(C1R“  C  x  (t|T)+C!R“  C  X 
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By  symmetry,  we  can  also  transform  the  layers  associated  to 

~(2)  I 

x  (t|T).  Then,  by  superposition,  we  see  in  Figure  4.2a  that  the 
internal  variables 


£(t)  -  x(t |t)  -  X(1) (t | T)  -  x(2)  £ t I T) 
5(t)  -  X(t|T)  -  X(1)(tjT)  -  X(2)(t|T) 


(4.6) 


can  be  obtained  by  composing  the  global  scattering  layers  with 
internal  sources 


q(t)A 


—  I  A  (  1  )  | 

'C2R2  °2X  (tlT) 


+  c^r‘1c1x(2) (t |t) 


(4.7) 


where  q(t)A  is  obtained  by  superposing  the  sources  of  the  global 
layers  with  the  ones  obtained  in  Figure  4.1c  for  the  local  layers.  The 
boundary  conditions  for  this  scattering  problem  are  also  obtained  by 
superposition  and  are  given  by 


C(0)  -  -x(0)  -  P(0)«5(0) 
6(T)  »  0 


(4.8) 


They  are  described  in  Figure  4 . 2b  where  we  have  represented  the 
aggregate  scattering  medium  associated  to  £(t)  and  5(t). 

To  solve  the  scattering  problem  satisfied  by  £(t),  we  can  use 
any  of  the  global  smoothing  algorithms  with  the  observations 
(yx(t) ,  y2(t))  replaced  by  -  (CjX(2) (t|T) ,  C2x(1)(t|T))  and  with 
initial  conditions  (-x(O) ,  P (0) ) .  The  two-filter  solution  of  this 


problem  is 
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+  P^CtiCgCt) 


(4.9) 


where 

5y(t)  -  (A-Pp(C^R*1C1+C^R“1C2))5F(t) 

"PF(CiRlL<:i*<2)  (tlT)+C2R2lc2*<1)  (tlT)) 
^(0)  -  -x(0) 


(4.10) 


and 

'PB(Ciai1VU>  <t|T)+C^R'1CJiUI  (t | T> )  (4<u 

^(tjCrW-o 

and  where  P  ,  P  and  P  denote  respectively  the  global  smoothing  and 

b  s  O 

forward  and  backward  filtering  error  covariances.  They  are  given  by 


P'V)  -  P^Jt)  +  p’X(t)  (4.12) 

with 

Pp(t)  -  APp+PpA'+BB'-PptC^R’^+C^R'^JPj,  (4.13) 

Pp(0)  »  P(0) 

and 

-p_ (t)  -  -ap -p_a,+bb,-p_(c'r71c..-k:;r"1c  )p  (4.14) 

p"1^)  -  0. 

0 

Since  the  equations  (4.9)- (4.11)  for  £(t)  depend  only  on  the  local 

*  (1)  i  *  (2)  i 

smoothed  estimates  x  (t|T)  and  x  (t|T5 ,  and  not  on  the  measurements 
y^  ( • )  and  y2 ( • ) ,  they  can  be  used  by  the  central  processor  to  compute  the 

centralized  estimate 

x(tjT)  -  5(t)  +  x(1)(t|T)  +  x(2)(t|T). 
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This  algorithm  was  first  derived  by  purely  algebraic  methods  in 
[11] .  One  interesting  application  of  this  result  is  for  the  map 
combining  problem  in  the  study  of  random  fields.  In  such  a  problem,  we 
consider  a  random  field  (gravitational,  magnetic,  or  else)  which  is 
modeled  by  the  process  x(t)  considered  above,  and  we  are  given  two  maps 
x^(t|T)  i*l,2  obtained  after  two  different  surveys  of  this  field. 
Then,  one  wants  to  combine  these  maps  to  obtain  a  global  map  x(t|T). 
This  cam  be  done  by  using  the  previous  algorithm  (see  [11]  for  more 
details) . 


B.  Decentralized  filtering 

To  obtain  some  similar  formulas  for  the  filtering  estimates  x(t|t), 
x^(t|t)  and  x*2*(t|t),  we  need  to  apply  the  superposition  principle  to 
the  transmission  layers  associated  to  the  filtering  problem  in  information 
filter  form.  Let  P(t)  and  (t)  i*l,2  be  the  solution  of  the  Riccati 
equations 


P  -  ap+pa'bb'-p^r’^+c^r"1^)? 


(4.15) 


and 


p(i)  *  ap^+p^a'+bb'-p^cirTVp^ 

XI  1 


(4.16) 


with  P(l)(0)  »  P(0)  i=l , 2 .  Then  if 


d(t) 


-1  ~  (i) 

P  (t)x(t)  and  dw  (t) 


P(iW  x(i)(t) 


(4.17) 


denote  the  global  and  local  filtered  estimates  of  x(t)  in  information 
form,  the  transmission  layers  associated  to  d  and  d^ 


acre 
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described  in  Figures  4.3a  and  b.  These  media  are  not  the  same,  but 
by  adding  and  subtracting  to  the  local  transmission  matrices,  we  can 
modify  the  local  layers  as  shown  in  Figure  4.3c,  so  that  the  trans¬ 
mission  media  are  now  identical.  By  superposition,  we  find  in  Figure 
4.4  that 

l(t)  -  d(t)  -  d(1)(t)  -  d(2)(t)  (4.18) 

satisfies  a  transmission  problem  with  the  same  medium  as  the  global 
filtering  problem  and  with  sources 

q(t)A  -  ((P(1,“1-p'1)BB'da)  +  {P(2)“1-p“1)B8,d(2,)A  .  (4.19) 

Thus,  l(t)  is  given  by 


l(t)  -  (A'+P^BB'JKt) 


+((P(1)“1-p"1)BB'd(1)+{P(2)~1-p‘1)BB'd(2) ) 


with  initial  conditions 


(4.20) 


1(0)  -  -d(0)  -  -p'1(0)*(0)  . 


(4.21) 


If  instead  of  computing  l(t),  we  want  to  compute 


r (t)  -  P(t)l(t)  -  x (t)  -  P(t) (P(1)“1x(1)+P(2)'1x<2)) 


(4.22) 


we  find  that 
r(t) 


(A-P  <ciRi1(:1+C2R2Ic2) }  1  (t) 
.(l)-l 


+  (PPl  -I)B8,P<1}”Ix<1) (t)+(PP(2)“1-I)BB,Pt2)“1xt2) 


(4.23) 

(t) 
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with  r(0)  *  -x(Q) ,  an  algorithm  which  is  the  same  as  the  one  obtained 
algebraically  by  Speyer  19]  and  Willsky  et.al.  [11]. 


Implementation 

The  previous  decentralized  filtering  algorithm  can  be  implemented 
in  two  ways,  depending  on  whether  1 (t)  or  r(t)  are  computed  in  a 
centralized  or  decentralized  way,  as  shown  in  Figures  4.5  and  4.6 
(see  also  Chong  [10]). 

In  the  centralized  method,  r(t)  is  computed  by  the  central  pro¬ 
cessor,  and  in  the  decentralized  method  the  local  processor  computes 


rU)  (t)  =»  (A-P(C|R~1C1+C^R"1C2))rU>  (t) 

+  (PP<l)~1-I)SB,P(l)"1x(l) (t) 

with  i~l,2.  The  initial  conditions  can  be  chosen  such  that 
r(1)  (0)  =  -x(0),  r(2)  (0)  -  0  • 


(4.24) 


(4.25) 


Then,  by  linearity,  one  gets 

r (t)  -  r(1) (t)  +  rl2) (t> 


(4.26) 


(i )  (i) 

so  that  if  each  local  processor  transmits  (x  ,r  ),  the  central 

^  (i) 

estimate  x  is  obtained  by  a  static  linear  combination  of  the  x  ' s  and 
r(i) 's,  i.e. 

x«r(1>  +r<2)  +p(p<1>-^>+P<2)-1;<2))  . 


(4.27) 
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(4.27) 


The  advantage  of  such  a  procedure  is  that  if  the  communications  between 
the  local  processor  and  the  central  one  are  restricted;  and  if  the 

/s 

central  processor  needs  to  compute  x  only  from  time  to  time,  then  the 

a  ( n  /  n 

local  processors  can  send  (x  , r  )  at  these  times  instead  of  sending 


The  advantage  of  such  a  procedure  is  that  if  the  communications  between 


the  local  processor  and  the  central  one  are  restricted;  and  if  the 

/s 

central  processor  needs  to  compute  x  only  from  time  to  time,  then  the 

a  ( n  Li.) 

local  processors  can  send  (x  ,r  )  at  these  times  instead  of  sending 
(i) 

x  continuously.  However,  a  disadvantage  of  this  method  is  that  it 

assumes  that  each  local  processor  knows  completely  the  global  model,  i.e. 

C.  and  R.  for  jj<i. 

2  2 


V.  The  smoothing  update  and  real-time  smoothing  problems 

In  some  cases,  instead  of  being  given  the  local  smoothed  estimates 

~ (i)  |  a  . 

x  ( 1 1 T }  i*»l,2  to  compute  x(t|T),  one  is  given  a  first  estimate 

Mi)  i 

x  ( t | T)  and  some  additional  measurements  y2(t),  0<t<T.  Then,  instead 

A  . 

of  reprocessing  all  the  data  (y^ < • ) ,y2 (• ) )  to  compute  x(t|T),  one  is 

^  (i)  | 

interested  in  updating  the  first  estimate  x  (t|T)  to  incorporate  the  new 
measurements  y2<-).  This  problem  will  be  called  in  the  following  the 
smoothing  update  problem.  It  arises  for  example  in  the  analysis  of  random 
fields  [11] ,  when  one  wants  to  update  an  old  map  to  incorporate  the  results 
of  new  surveys. 

A  variant  of  this  problem  that  will  also  be  discussed  in  the  following 
is  the  real-time  smoothing  problem.  In  this  case,  if 

Yx  *  H(m;y^(s) ,  0<s<T) 

Y  2  “  H(m;y  (•),  0<s<t) ,  0<t<T 
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denote  the  Hilbert  spaces  spanned  by  all  the  observations  ( • )  and 
by  the  observations  y2 ( * )  up  to  time  t,  one  wants  to  compute 

^(t)  -  Etx(t)  (5.1) 

A  , 

in  terms  of  x  (t|T)  and 


A.  The  smoothing  update  problem 

By  superposing  the  global  and  local  scattering  layers  obtained 
for  the  decentralized  smoothing  problem,  and  by  denoting 


H(t)  -  x(t|T)  -  i(1)(t|T)  (5.2) 

8(t)  -  X(t|T)  -  X(1)(t|T)  (5.3) 


we  obtain  the  scattering  layers  described  in  Figure  5.1.  This  shows 
that  the  variables  n(t)  and  9 (t)  can  be  obtained  by  composing  the 
global  scattering  layers  with  internal  sources 


q  (t)  A 


( 

'‘•CSR2l(VV  11 


A 


(5.4) 


The  boundary  conditions  for  this  scattering  problem  can  also  be 
obtain  by  superposition,  and  they  are  given  by 

n(0)  -  -P (0)9(0) ,  8(T)-0  .  (5.5) 

To  solve  this  scattering  problem,  we  can  use  again  the  two- filter 
formula.  We  obtain 

p’^tjntt)  -  p"1(t)np(t)  +  p^fongft) 


(5.6) 
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Figure  5.1.  Scattering  layers  associated  to  the  smoothing 
update  problem. 


Figure  5.2a.  Scattering  medium  for  the  real-time 
smoothing  problem. 
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Figure  5.2a.  Scattering  medium  for  the  real-time 
smoothing  problem. 
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where 

ny(t)  -  (A-pp(c^R"1c1+c^R“1c2))np(t) 

+  PpC^R"1(y2(t)-C2x(1) (t|T)) 

1^(0) -0 

and 


+  PBC2R21(y2(t)“C2*(1) (tlT>) 

^  (T)  TTfl  (T)»0  , 


and  where  P  .  P  and  P  are  given  by  (4.12)- (4.14) . 
S  F  3 


(5.7) 


(5.8) 


B.  The  real-time  smoothing  problem 

The  aggregate  scattering  medium  associated  to  the  real-time 

smoothing  problem  is  given  by  Figure  5.2a  where  the  scattering  matrix 

S(t,0)  and  the  sources  q(t,0)  are  generated  by  the  global  scattering 

layers  described  in  Figure  4.1a  and  where  the  medium  S  ^  (  t, t)  and  the 

sources  q  ^  (T,t)  are  generated  by  the  local  layers  of  Figure  4.1b. 

-'(I)  i 

Similarly,  the  local  scattering  medium  associated  to  x  (t|T) 
is  given  by  Figure  5.2b,  where  the  right  layer  is  the  same  as  for 

A 

x  (t) .  The  medium  of  the  left  layer  cam  be  made  the  same  as  for 
x  (t)  if  we  use  the  modification  of  the  local  layers  described  in 
Figure  4.1c.  By  doing  so,  we  find  that 


-38- 


where 

ny(t)  -  (A-pp(c^R"1c1+c^R“1c2))np(t) 

+  PpC^R"1(y2(t)-C2x(1) (t|T)) 

1^(0) -0 

and 


+  PBC2R21(y2(t)“C2*(1) (tlT>) 

^  (T)  TTfl  (T)»0  , 


and  where  P  .  P  and  P  are  given  by  (4.12)- (4.14) . 
S  F  3 


(5.7) 


(5.8) 


B.  The  real-time  smoothing  problem 

The  aggregate  scattering  medium  associated  to  the  real-time 

smoothing  problem  is  given  by  Figure  5.2a  where  the  scattering  matrix 

S(t,0)  and  the  sources  q(t,0)  are  generated  by  the  global  scattering 

layers  described  in  Figure  4.1a  and  where  the  medium  S  ^  (  t, t)  and  the 

sources  q  ^  (T,t)  are  generated  by  the  local  layers  of  Figure  4.1b. 

-'(I)  i 

Similarly,  the  local  scattering  medium  associated  to  x  (t|T) 
is  given  by  Figure  5.2b,  where  the  right  layer  is  the  same  as  for 

A 

x  (t) .  The  medium  of  the  left  layer  cam  be  made  the  same  as  for 
x  (t)  if  we  use  the  modification  of  the  local  layers  described  in 
Figure  4.1c.  By  doing  so,  we  find  that 


-41- 


and 

P“*(t)  =  P^tt)  +  P'1,_1(t)  ‘ 

RS  F  B 

From  (5.13)  one  get  ni^(t)=Of  so  that  (5.12)  reduces  to 

3 

*RS{t)  =  *(1)  (t i t)  +  . 


(5.15) 


(5.16) 


A  similar  problem  to  those  discussed  here  appears  when  one 
considers  a  decentralized  filtering  structure  where  instead  of  assuming 
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and 


(6.4) 


(6.5a) 


(6.5b) 


A  I 

With  this  model  structure,  the  smoothed  estimates  x(t|T)  and 

*  (i)  . 

(t|T)  can  be  obtained  by  composing  the  global  and  local  scattering 
layers  described  in  Figures  6.1a  and  6.1b.  In  Figure  6.1a,  it  is 
assumed  that  the  matrices  AjB,^  and  C 2  have  the  same  block  structure 
as  in  (6. 4)- (6. 5).  The  main  feature  of  the  local  scattering  layers  is 
that  the  internal  variables  x^  (t J  T)  and  A^  (t  |  T)  are  of  smaller 
dimension  than  the  variables 

/  xL{t|T>  \  /x^tltA 

x(t|T)  =  (  )  ,  A(t|T)  4  )  (6.6) 

\s1(tjT)/  \  CTX (t|T)  / 


associated  to  the  global  scattering  problem.  However,  by  adding  some 
appropriate  internal  sources,  the  local  scattering  layers  can  be 
imbedded  in  the  global  ones  as  shown  in  Figure  6.1c  and  the  internal 


variables 
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simplified  models  which  describe  only  the  evolution  of  a  set  of 
local  variables  (machine  angles  and  frequencies) .  This  motivates  us 
in  this  section  to  consider  the  case  when  the  global  model  is  given  by 
(4.1)  and  (4.2)  and  when  the  local  models  and  measurements  are  given 
by 

\  am  21  v  /f  \  X  Q  n 


*21^  (t|T)+B21B^X^1^  (t|  T) )  A 


0 
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variables 
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/  (t|T)\ 

^  \  °  / 

have  now  the  same  dimension  as  x(tjT)  and  X(tjT). 

By  superposing  the  global  scattering  layers  of  Figure  6.1a  and 
the  imbedded  local  layers  of  Figure  6.1c,  we  obtain  the  scattering 


*21^  (t|T)+B21B^X^1^  (t|  T) )  A 


0 
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seek  below  to  transform  the  problem  so  that  B^B^-O.  When  this  is 
the  case,  the  two- filter  solution  of  the  scattering  problem  is  given 
by 


and 


<-a-pb (c;R-11c1«’r-1c2  )) 


5b") 


nB(t) 


(6.11) 


M> 
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( 


and 


have  now  the  same  dimension  as  x(tjT)  and  X(t|T). 

By  superposing  the  global  scattering  layers  of  Figure  6.1a  and 
the  imbedded  local  layers  of  Figure  6.1c,  we  obtain  the  scattering 
layers  described  in  Figure  6.2.  This  shows  that  the  internal  variables 


and 


S(t) 


'^(tjT)  -  (tjT) 


n(t) 


i^tlT) 


X1(t|T)  -  X^1} 
O^tjT) 


(6.7a) 


(6.7b) 


satisfy  a  scattering  problem  with  the  same  medium  as  for  the  global 
smoothing  problem  and  with  internal  sources 


and 


3+<t  )A 

q_(t)A 


0  V 

(t|T)  -  B21B[^L)  (tjT)  J 
R21(Y2(t)  “  C21*l1} (tlT))A  ’ 


(6.8a) 


(6.8b) 


These  sources  require  the  knowledge  of  X(^  (t|T)  as  well  as  of 
(tjT) .  This  is  not  very  convenient,  and  consequently  we  will 
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and  where  P  ,  P  and  P  are  given  by  (4. 12}- (4. 14)  with  the  special 
S  F  B 

matrix  structure  (6.4)-(6.5).  Since  equations  (6. 10)- (6. 11)  depend 

-  (1)  , 

only  on  the  local  estimates  x^  (t|T)  and  on  the  new  measurements 
y2  ( • ) ,  we  can  use  them  to  compute  the  centralized  smoothed  estimate 


To  show  that  we  can  always  transform  the  problem  so  that 
B2jBj=0,  we  consider  the  process 

^u(t)  .  (6.13) 

By  projecting  w2(t)  on  w^t),  one  gets 


w2(t)  »  E [w2 (t)  iw^ (t) ]  +  w2(t) 


TlWi(t,  +  w2(t) 


so  that 


0  -  Elw  (t)w^(t)]  *  E[w2<t)w|(t)]  -  T^E [w^ (t)w^(t)] 

*  B  B'  -  T  B  B*  . 

*211  ill 


Then,  if  we  consider  the  similarity  transform 


CM; :)(::) 


(6.14) 


(6.15) 


(6.16) 


C21R2  (VC21X1  (t'T))A 


Scattering  layers  obtained  by  superposition 


-50- 


the  system  (6. 4) -{6. 5)  becomes 


where  we  have  B^B*'  »  (B2l“,riBl*Bi*0,  so  that  we  can  apply  the 
formula  (6- 9) -6. 11)  to  this  system. 

Some  similar  results  cam  be  obtained  for  the  filtering  estimates. 

The  problem  that  we  shall  study  is  one  where  the  processor  can  com¬ 
municate  with  S2  and  where  the  model  available  to  S2  is  the  global  model, 
whereas  the  model  available  to  is  of  smaller  size.  In  this  case,  the 

a  (2) 

estimate  (t)  (based  on  the  observations  y^f*))  is  received  by  the 

processor  S2  which  then  wants  to  combine  this  estimate  with  the 

A 

observations  y2  M  in  order  to  obtain  x(t).  This  estimation  structure  is 
not  synnetric,  but  in  this  case  the  decentralized  filter  has  a  simple  form 
which  would  otherwise  be  hidden  by  several  changes  of  basis  for  the 


global  state-space. 
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seek  below  to  transform  the  problem  so  that  B^B^-O.  When  this  is 
the  case,  the  two- filter  solution  of  the  scattering  problem  is  given 
by 


and 


<-a-pb (c;R-11c1«’r-1c2  )) 


5b") 


nB(t) 


(6.11) 
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tn  the  following,  we  will  denote  respectively  by  P(t)  and  (t) 

the  solutions  of  the  Riccati  equations  (4.15)  and 


,111  .  a  p'1'  ♦  p^'a'  ♦  s  a;  -  P^VjVI11  . 


(6.18) 


where  in  (4.15)  we  assume  that  A,B,C^  and  have  the  block  structure 
(6. 4) -(6. 5).  The  initial  conditions  for  these  equations  cure  such  that 


if 


we  have  P^(0) 


Then  if  we  denote  by 


d(t) 


(6.19) 


(6.20) 


the  information  filter  form  of  the  centralized  filtering  estimate,  the 

global  transmission  layers  associated  to  d(t)  are  described  in 

Figure  6.3a.  Similarly,  if  d^1^  (t)  =  1  (t)x^  (t)  ,  the  local 

transmission  layers  associated  to  d|^  (t)  are  given  by  Figure  6.3b. 

Since  d^  (t)  is  of  smaller  size  than  d(t),  we  cannot  apply  directly 

the  principle  of  superposition  to  the  global  and  local  layers.  However, 

as  for  the  smoothing  problem,  we  can  imbed  the  local  transmission  layers 

(  1  )  *  * 

inside  the  global  ones  as  shown  in  Figure  6.3c.  The  vector  (d^  (t),0) 

of  transmitted  variables  has  now  the  same  size  as  d(t).  Then,  by 


i 

i 
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,HiRIlyi*CilB2ly2>a 


0— — ►  <vt+A> 


(t+A) 


C22R2  y2A 


i  -  (A^p^1,"1B1a^)A 


d^1) (t+A) 


HiRIlyia 


B  B' 

12  1 


(a)  Global  transmission  layers,  (b)  Local  transmission  layers, 
(c)  Imbedding  of  the  local  layers. 
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superposing  the  transmission  layers  of  Figures  6.3a  and  6.3c,  we  find 
in  Figure  6.4  that 


l(t)  =  d(t) 


(6.21) 


satisfies  a  transmission  problem  with  the  same  medium  as  for  the 
global  filtering  problem  and  with  sources 


q(t)A 


(t)A  +  C^R'S^A 


(6.22) 


(note  that  b12bi  =  TiBiBi^ ‘  Thus,  (t)  is  given  by 


l(t)  =  - (A'+P~  bB‘ ) l(t)  +  C'b"1 Y2(t) 


with  initial  conditions 

/pj1)'1(0)x1O) 
1(0)  =  P~1(0)x(0)  -  l 

\  0 


If  instead  of  computing  l(t)  we  want  to  compute 


r(t)  =  P(t)  l(t)  *  x(t) 


(6.23) 


(6.24) 


(6.25) 


0 
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superposing  the  transmission  layers  of  Figures  6.3a  and  6.3c,  we  find 


in  Figure  6.4  that 


l(t)  =  d(t) 


(6.21) 


satisfies  a  transmission  problem  with  the  same  medium  as  for  the 


without  revealing  their  position.  These  units  tend  consequently  to 
communicate  as  rarely  as  possible,  in  addition,  seme  delays  are  incurred 
when  transmitted  signals  need  to  be  relayed  through  several  communication 
nodes  before  reaching  a  processing  center.  Some  noise  appears  also  in 
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without  revealing  their  position.  These  units  tend  consequently  to 
communicate  as  rarely  as  possible,  in  addition,  seme  delays  are  incurred 
when  transmitted  signals  need  to  be  relayed  through  several  communication 
nodes  before  reaching  a  processing  center.  Some  noise  appears  also  in 
communications  when  the  environment  interacts  with  the  transmitted  infor¬ 
mation  through  distortions,  background  noise,  multipath  effects, 
jamming,  etc. . 

The  extension  of  the  results  of  the  previous  sections  to  take  into 
account  the  existence  of  delays  and  of  noise  in  the  communication  channels 
is  presently  under  consideration. 

We  note  also  that  in  this  paper,  no  special  assumption  has  been  made 
on  the  structure  of  the  dynamical  systems  considered.  When  there  exists 
a  natural  separation  of  the  local  systems  arising  either  from  a  weak 
coupling  of  these  systems,  or  from  a  separation  between  fast  and  slow  time 
scales,  we  can  use  the  scattering  framework  to  obtain  same  decentralized 
estimation  algorithms  which  exploit  this  structure,  as  will  be  shown  in  a 
subsequent  paper. 
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